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We study the adsorption of homogeneous or heterogeneous polymers onto heterogeneous planar 
surfaces with exponentially decaying site-site correlations, using a variational reference system ap- 
proach. As a main result, we derive simple equations for the adsorption-desorption transition line. 
We show that the adsorption threshold is the same for systems with quenched and annealed disorder. 
The results are discussed with respect to their implications for the physics of molecular recognition. 



I. INTRODUCTION 

Polymer adsorption has been studied intensively in the 
last decades because of its wide practical applications and 
rich physicsiiS. The early studies mostly focused on the 
adsorption onto physically and chemically homogeneous 
surfaces. In nature, however, real objects are very rarely 
perfectly homogeneous with ideal geometry. Real sub- 
strates usually contain some heterogeneity. Therefore, a 
number of more recent studies have considered the influ- 
ence of random heterogeneities on the adsorption transi- 

From a more fundamental point of view, the adsorption 
of random heteropolymers (RHPs) onto heterogeneous 
surfaces has attracted much attention, because RHPs 
serve as biomimetic polymers to study the physics of pro- 
teins and nucleic acids. Studying the interaction of RHPs 
with random heterogeneous surfaces can provide insight 
into the physics of molecular recognition, a process play- 
ing a key role in living organisms (e. g., enzyme-substrate 
or antigen- antibody interactions). 

The present work deals with the adsorption of ho- 
mopolymers and heteropolymers on surfaces with chem- 
ical heterogeneities. As a general rule, the introduction 
of chemical heterogeneities at constant average load pro- 
motes adsorption, since the polymers can adjust their 
conformations such that they touch predominantly the 
most attractive surface sites. The adsorption transition is 
thus shifted to a lower value of the mean polymer/surface 
affinity. The extent of the shift is determined by the char- 
acteristics of the randomness. In the present work, we 
focus on the influence of correlations on the adsorption 
transition. 

Our study was motivated by the idea that clustering 
and the interplay of cluster sizes should contribute to 
molecular recognition. This idea is supported, among 
other, by a series of recent papers by Chakraborty 
and coworkera^iiiSiSiLSiSiiSiii. Using different analyt- 
ical approaches and Monte Carlo (MC) simulations, 
these authors have established that below the adsorption 
threshold, statistically blocky polymers favor statistically 
blocky surfaces over uncorrelated surfaces. They found 
that, upon increasing the strength of the interactions, 
the adsorption transition of heteropolymers on heteroge- 
neous surfaces is followed by a second sharp transition, 
where the polymers freeze into conformations that match 
the surface oattern. This haDDens at adsorntion energies 



signiflcantly beyond the adsorption transition. In the bi- 
ological context, a two-step adsorption process where ad- 
sorption is followed by "freezing" possibly describes the 
physics of protein/DNA recognition, where the protein 
slides on the DNA for a while before finding it's specific 
docking site. For other molecular recognition processes, 
such as antibody-antigen interactions, it is not only im- 
portant that a protein adsorbs to a given particle, but 
equally crucial that it does not adsorb to other objects. 
In order to assess the importance of cluster size matching 
in those cases, one must calculate the shift of the adsorp- 
tion transition as a function of correlation lengths on the 
polymer and on the surface. 

In a previous publication^, we have calculated the 
adsorption/dcsorption transition of an ideal (non-self- 
interacting) RHP with sequence correlations on an im- 
penetrable planar homogeneous surface. We have com- 
pared two different approaches. First, we have calculated 
numerically the partition function of a random sample of 
lattice RHPs^. Provided the sample is sufficiently large, 
this approach is basically exact. Second, we have applied 
a reference system approach which was originally devel- 
oped by ChenS^ and Denesyuk and Erukhimovich^S in 
a similar context (RHP localization at liquid-liquid in- 
terfaces). The second approach allowed to derive an ex- 
plicit and surprisingly simple expression for the adsorp- 
tion transition point (repeated here in Eq. (|42|) '). which 
fitted the numerical results excellently. 

Having thus established the validity of the reference 
system approach, we use it here to study the adsorption 
of a single Gaussian chain onto a heterogeneous planar 
surface with exponentially decaying site-site correlations. 
We consider both homopolymers and random heteropoly- 
mers. As in the case of the RHP on an homogeneous 
surface, we shall derive simple analytic formulae for the 
adsorption/dcsorption transition point. This will allow 
to discuss how surface patchyncss infiuences the adsorp- 
tion behavior, and whether the interplay of surface and 
polymer correlations can bring about something relevant 
to molecular recognition in the vicinity of the adsorption 
transition. 

The problem of adsorption on heterogeneous surfaces 
has already been considered theoretically by a number 
of authors, both for homopolymers and heteropolymers 
(block-copolymers, periodic or random copolymers). 

Baumgartner and Muthukumai— studied the effect of 
chemical and physical surface roughness on adsorbed ho- 



2 



mopolymers. Using MC techniques, they calculated the 
adsorbed homopolymer amount on chemically heteroge- 
neous surfaces, as a function of temperature and deter- 
mined the shift of the adsorption-desorption transition 
temperature with decreasing concentration of adsorbing 
sites for different chain lengths. In addition, they pre- 
sented a scaling analysis showing that lateral size of the 
adsorbed chain depends on the distribution of re- 
pelling "impurities" on the adsorbing surface and the 
strength of their repulsive interaction with the polymer. 

is unaffected by randomness and simply equal to the 
radius of gyration of a Gaussian chain in two dimensions 
if impurities are very rare or weak. However, when the 
impurities are not rare and repulsive barriers are high, 
i?l I is of order of the mean distance between neighboring 
impurities. 

Similar conclusions were reached by Sebastian and 
Sumithraiiii, who developed an analytical theory of the 
adsorption of Gaussian chains on random surfaces, using 
path integral methods combined with the replica trick 
and a Gaussian variational approach. They took surface 
heterogeneity into account by modifying de Gennes' ad- 
sorption boundary condition^Si^, and analyzed the influ- 
ence of randomness on the conformation of the adsorbed 
chains. 

MC simulations of single polymer adsorption have also 
been performed by Moghaddam and WittingtoniSi. The 
sites on the heterogeneous surface were taken to be sticky 
or neutral for the polymer. The temperature dependen- 
cies of the energy and the heat capacity were calculated 
for different sticker fractions and compared for the cases 
of quenched and annealed surface disorder. Huber and 
Vilgisi^ have studied the case of physically rough sur- 
faces. 

The above mentioned works dealt with the single-chain 
problem. A number of other authors investigated the 
adsorption from solution onto heterogeneous surfaces us- 
ing MC simulationsiLi^, numerical lattice self-consistent 
field (SCF) calculationaiS, and an analytical SCF thc- 
orjiSS. 

As mentioned earlier, the problem of heteropolymcr 
adsorption on heterogeneous surfaces has been studied 
extensively by A. Chakraborty and coworkers. Most of 
the results were summarized in a recent review articleflTl 
Srebnik et al. used replica mean field calculations and the 
replica symmetry breaking scheme to study the adsorp- 
tion of Gaussian RHP chains for theta solvents^ and for 
good solventsi^. MC simulations combined with simple 
non-replica calculations allowed to analyze the role of the 
statistics characterizing the sequence and surface site dis- 
tribution, and that of intersegment interactions, for the 
pattern recognition between a disordered RHP and a dis- 
ordered surface^i^iSiL2iS. 

The adsorption of globular RHPs (e. g., RHPs in a bad 
solvent) on patterned surfaces has recently been consid- 
ered by Lee and Vilgis with a variational approach^i. A 
MC simulation study of single RHP chain adsorption onto 
random surfaces has also been performed by Moghad- 
damS^. 



Interesting applications of interactions between copoly- 
mers and heterogeneous surfaces have been examined 
by Genzeri2i2li2Si2&. Using a three dimensional self- 
consistent field model, he studied the adsorption of A- 
B copolymers from a mixture of A-B copolymers and 
A homopolymcrs onto a heterogeneous planar substrate 
composed of two chemically distinct regularly or ran- 
domly distributed sites. He demonstrated that via sub- 
strate recognition, the copolymer can transcribe the two- 
dimensional surface motif into three dimensions. The 
specific way the surface motif is translated was shown 
to be strongly dictated by the copolymer sequence. 

The present work differs from all these approaches in 
that we focus on the adsorption transition, i. e., on the 
region in parameter space where adsorption just sets in, 
and that we derive simple analytical expressions for the 
location of that transition. We shall see that the char- 
acteristic features observed in the regime of stronger ad- 
sorption do not necessarily carry over to this regime of 
weak adsorption. 

The rest of the paper is organized as follows: Sec. ^ 
introduces the model. The free energy calculation based 
on the reference system approach is explained and carried 
out in Sec. IIIII Some aspects of the calculations in this 
section are very similar to those already described in our 
earlier publication, Ref. [23. Therefore, we will skip the 
corresponding intermediate technical steps and refer the 
reader to Ref. for more details. Explicit analytic ex- 
pressions for the adsorption transition point are obtained 
in Sec. IIII El We summarize and discuss our final ex- 
pressions in Sec. IIVI Some intermediate calculations are 
described in three Appendices. 



II. THE MODEL 

We consider single Gaussian polymers chain consisting 
of N monomer units near an impenetrable planar surface. 
The surface carries randomly distributed sites of differ- 
ent chemical nature. The surface pattern is described 
by a locally varying interaction parameter cr(x). Here x 
denotes Cartesian coordinates {x,y). A particular (het- 
ero)polymer realization is described by a sequence ^(n), 
where n indicates the monomer number [n € [0 : N]). 
The strength and the sign of the monomer-surface inter- 
action is determined by the product £_{n)a{x.) (a "-f-" sign 
refers to attraction, "-" to repulsion). 

The Hamiltonian of the system can, therefore, be rep- 
resented as follows: 

kf-2^1 '-[o-n] ^-^^(-M-MxW), 

(1) 

where r(n) is the chain trajectory in the space, r = {x, z}, 
(3=1/ kT is the Boltzmann factor, a is the Kuhn segment 
length of the chain, and V{z) describes the shape of the 
monomer-surface interaction potential. This potential is 
taken to be a delta-function pscudopotential shifted at a 
small but finite distance zq from the impenetrable surface 
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and has the form 



V{z) 



6{z - zq) if > 
+00 if z < 



(2) 



The model Hamiltonian |^ can be used to describe 
interactions of both heteropolymers and homopolymers 
with both heterogeneous and homogeneous surfaces. Ho- 
mopolymers are described by a constant sequence func- 
tion, i.e. ^(n) = ^0- The sequences ^(n) of heteropoly- 
mers are taken to be randomly distributed according to 
a Gaussian distribution function 



P{^{n)} cx 



(3) 



exp 



N 



dni o^7^2(^(?^l)-Co) c (^1,712) iC{n2)-^o) 



where ^0 = (0 characterizes the "mean charge" of the 
chain, c(ni,n2) describes sequence correlations, and 
is the inverse of c, defined by Jdnc~'^{ni,n)c{n,n2) = 
d{ni — 712). Here we consider exponentially decaying cor- 
relations, i. e.. 



c(ni,n2) 



mni) 

c{\n2 - 



-eo)(f(n2) 



-eo)) 

-rp|?i2- 



(4) 



rill 



The parameter Tp corresponds to an inverse "contour cor- 
relation length" on the polymer, and Ap gives the vari- 
ance of the single-monomer distribution. The propor- 
tionality factor in Eq. (|3Jl is chosen such that P{£^{ti)} 
is normalized. This heteropolymer model also covers ho- 



where (Tq ~ {<^) is the mean load of the surface, and the 
site-site correlation function g{x.i,x.2) decays exponen- 
tially: 

g(xi,X2) = ((cr(xi) - cro)(CT(x2) - ctq)) (6) 
= .9(1x1 -X2|) = A2e-r=l'^i-''^l. 

Here is again the variance of the single-site distri- 
bution on the surface and is the inverse correlation 
length of the surface pattern. As before, g~^(xi,X2) 
is the inverse function of (7(xi,X2) with the property 
/ dx (7~^(xi, x)g(x, X2) ~ (5(xi — X2), and P{(t(x)} is nor- 
malized. Homogeneous surfaces are obtained in the limit 
As 0, and are characterized by cr(x) = ctq- 



III. THE REFERENCE SYSTEM APPROACH 

A. Quenched Average via Replica Trick 

The free energy of the system is given by the logarithm 
of the conformational statistical sum averaged over all 
realizations of the disorder f3F = (hiZ){o.(x),^(n)}- Here 
and below angular brackets (. . .) stand for disorder aver- 
age. To perform the quenched average we use the well- 
known replica trick 



mopolymers in the limit 



0. 



Similarly, the distribution of surface sites on heteroge- 
neous surfaces is chosen to be Gaussian, with the proba- 
bility 



(InZ) 



{<T(x),C(n)} 



= lim 

?n — >0 



1 



{<t(x),5(»)} 



(7) 



P{<7(X)} « 



(5) 



exp -| //dxi dx2 (cr(xi)-(To) 5 ^(xi,X2) (cr(x2)-CTo) 



First we average over the surface disorder o'(x). Intro- 
ducing 771 replicas of polymer chains and making use of 
the properties of the Gaussian distribution, Eq. (0), we 
obtain 



3 " fdr^ 



.-1 "'O 



/IIL 
l[ V{r^{n)} exp 

/3(To / dn£,{n)y^V{za{n)) + - II dni dn2 £.{ni) ^{n2)g{ni,n2) 
Jo ^ JJq 



(8) 



X exp 



where a is the replica index and / 'D{r{n)} denotes inte- 
gration over all possible trajectories r{n) of the chain. In 
Eq. we have introduced the notation 

?(77l,772)= (9) 
m 



for clarity and future reference. 

Next we average over the sequence distribution S,{n). 
This affects only the second exponential in ©. Again, 
exploiting the properties of Gaussian integrals, we obtain, 

m 
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){a(x),^(n)} 



y l[V{r^{n)} exp 



X cxp 



X cxp 



X cxp 



3 ^ /9r„ 



1 



Vdctl? 



(10) 



TV 



N 



CTo^o // dnidn2K ^{ni,n2) (3^V{za{n2)) 
dni dn2 dns K"^{ni,n2) g(ri2, "s) 

dnidn2dn3K^^{ni,n2)c{n3~ni) P"^ ^ ^(^^(ria)) ^(2^(713)) 



a,7— 1 



r 



where if (711,712) is defined as 

K{ni,n2) = 5{ni-n2)- I ^7733(711, 713)0(713-712) . (11) 
Jo 

Eq. p(J|) covers the cases of homogeneous and hetero- 
geneous polymers adsorbing on homogeneous and het- 
erogeneous surfaces. If at least one of the partners is 
homogeneous, the matrix K becomes unity, if (711,712) = 
5(711 — 712), anddetif = 1. If the polymer is homogeneous 
(a homopolymer), the argument of the last exponential 
of Hl()|l vanishes and we recover Eq. JHJ in the special case 
^(71) = ^0- If the surface is homogeneous, the argument 
of the second last exponential in (|10|l vanishes and we 
recover Eq. (6) of Ref. 

So far the discussion has been fairly general. In the 
following, we shall mostly restrict ourselves to two spe- 
cific situations: Homopolymer adsorption on heteroge- 
neous surfaces (c(7i) = 0, A'(7ii, rt2) = i5(7ii — 712)), and 
the adsorption of overall neutral heteropolymers on over- 
all neutral heterogeneous surfaces (ctq = Co = 0). In the 
second case, the last three exponentials in Eq. (|10|) can 
be dropped, and the effect of the disorder comes in solely 
through the determinant det K ^1. 



B. Definition of the Reference System 

To proceed with the calculations, we introduce a refer- 
ence systemS^ with conformational and thermodynamic 
properties reasonably close to those of the original sys- 
tem. A natural choice for the reference system is a sin- 
gle homopolymer near an adsorbing homogeneous sur- 
face, interacting via an attractive potential of the same 
form V{z) as in |j2J), with the interaction strength wq = 
"■q ^0 + 7^0- The parameter uq is a variational parameter 
which can be adjusted such that the reference system is 
as close as possible to the original one. The Hamiltonian 
of the reference system has the following form: 



kT 



3 
2^ 



on 



N 



l3wo / dnV{z{n)). (12) 



The next steps are very similar to those already de- 



eropolymcr chain adsorption on a homogeneous planar 
surface, and shall be only sketched very briefly here. We 
begin with recasting the quantity (Z™) from Eq. 1)1 Of) in 
the form 



(13) 



where Zq is the partition function of the reference system. 



and 



r i(m) 



denotes the chain average with respect to 771 



independent reference chains. The general expression for 
A,n is complicated, but can be obtained from Eq. (|10|l 
in a straightforward way. In the case of homopolymers 
adsorbed on a heterogeneous surface, the quantity Am is 
given by 



£0 
2 



N pN 

dnidn2^{ni , 712)— /3f 



dn^V{za{n)) . 

a=l 



(14) 

and in the case of a neutral heteropolymer on a neutral 
surface, we get 



A.„, = - - ln(det K) - (3vo / dii ^ V{z^{n)) 

•'^ a = l 



(15) 



The free energy difference between the original system 
and the reference system is given by: 



13 AF = lim ((Z") - Z")/77i 

lim([exp(A„)][,™)-l)/7n. 



(16) 



m-»0 



Unfortunately, the expression ( [exp (A„i )][,'"■' — 1) cannot 
be calculated exactly. Therefore, we expand it up to lead- 
ing order in powers of fg, g, and c (i. e., and Ap). In 
the homopolymer case, Eq. 1)14(1 . this is equivalent to ap- 
proximating exp(A) — 1 K. A. In the heteropolymer case, 
Eq. ((T3|l ') we also have to expand ln(det K): 

\n{detK) = Tr(ln A') w Tr(A' - 1) (17) 
dni dn2 5(111, ri2) c{n2 - 7ii) 
The resulting expression for homopolymers on a het- 
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/3AF 
N 



■pvo[Viz{n))]o (18) 



N 



+ fi^il j dk j dx' dx" 5(x' - x") [y(z(7i))(5(x(7i) - :&')V{z{n + fc))(5(x(n + fc) - x")]o 
- (3^ Co / t^x' dx" 5(x' - x") [y(z(n))J(x(n)-x')]o [V{z{n))5{-^{n) - x")]o ■ 

I 

The last term corresponds to contributions from different 

replicas, whereas the other two belong to single repli- r;^i ^ / X>{r(n)} cxp[-/3go(r(n))]X(r(n)) 

cas. Wc have used the fact that averages for different / P{r(n)} exp[— /3i/o(i"(fi))] 

replicas factorize, since different replicated chains are 

independent in the reference system. Thus [~]o denotes Likewise, we obtain for neutral heteropolymers on 
the average with respect to a single reference chain, neutral heterogeneous surfaces 



/3AF 
N 



Pvo[V{z{n))]o 



(20) 



N 



+ f3^ J dk c{k) J dx dx" g{x - x") [V{z{n))S{x{n) - x')V{z{n + k))S{x{n + k) - x")]o 
-P^ J dx dx" g{x - x") [V{z{n))S{x{n) ^ x')]o [y(z(n))(5(x(n) - x")]o ^ dk c{k) 



r 



As in Eq. H18|l . the first two terms contain contributions 
from the same replica, and the last term accounts for 
contributions from different replicas. 

At this point, the parameter vq is still a free variational 
parameter of the theory. In the final step, we will choose 
it such that the free energies of the reference system and 
the original system are the same, i.e. AF = 0. 



C. Green's Function and Transition Point of the 
Reference System 

To calculate the averages in (|18|l and (|20|l we need to 
find the Green's function of the reference system. It sat- 
isfies the following differential equation^ii^^: 

an 6 

with the appropriate initial and boundary conditions: 



G(r,r';0) = 6{r - r') , 
G(r,r';iV)||rHo, = 0, 
G(r,r';7V)U=o = 0. 



(22) 
(23) 
(24) 



First, one can separate the variables x, y and z writing: 
Giv,v';n)=GJx,^-n)GJz,z':n), (25) 



Gx(x, x'; 77,) is the free chain Gaussian distribution in two 
dimensions: 



Gx(x,x';77,) 



2'i:na 



exp 



3(x-x') 



l\2 



2na? 



(26) 



For the z-dcpendent part of the Green's function H25|l in 
the long-chain limit 77, ^ 1, one obtains in the ground- 
state dominance approximation^Si^ 



, fco/2 

l + (2fcb-M^o 



g-feblz'-zol _ g- 



(27) 

-ki(z' +zo) 



where kb, determining the ground-state energy (eo 
ki,a^/6), is the root of the transcendental equation 



2k 



I — g-2fczo 



ko with 



(28) 



The complete Green's function is calculated and dis- 
cussed in the Appendix IbI 

The adsorption transition in the reference system takes 
place when kb and eo vanish. This corresponds to the 
conditions^ zq^o = 1 or, according to the definition (|28|l 
for ko, to QBwo/a^ = \Izq. 
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D. Free energy difference between the reference 
system and the original system 

With the rcsuhs from the previous subsection, we 
can now evaluate the averages in (|18|) and (|20|l . Using 
the ground-state dominance approximation H27|) for the 
Green's function, we obtain the same resuh for [y(z(n))]o 
as in Ref. El 



[V(zin))]o = 



1 + {2kb - ko)zo 



(29) 



Correspondingly, the averages [V{z{n))S{:x.{n) — x')]o ap- 
pearing in the last term of (|18|l and (|20|) are given by 



[y(^(n))<5(x(n)-x')], 



2fc2/fco 



1 + (2fcf, - ko)zo 



(30) 



where S is the surface area. The additional factor 1/5 
in H3U|) . compared to H29|l . results from the translational 
freedom of the polymer in the surface plane. Substitut- 
ing this result into the last term of 1)18(1 and ((2()|l and 
performing the integration over / dx.' dx" gives the factor 
1/S. Since 5' oo, this contribution is infinitesimally 
small and can be omitted in the following calculations. 

This seemingly "technical" point has an important 
physical interpretation. After dropping the last term. 



Eqs. 1(18(1 and ((20(1 do no longer contain terms with con- 
tributions from different replicas. These terms, however, 
are the only ones which distinguish between quenched and 
annealed disorder. One can see that by directly carrying 
out an analogous calculation for annealed disorder. In 
the annealed case, one has to average directly over the 
partition function and calculate 



= ln(Z) -InZo =ln[exp(^i)]o. 



(31) 



After expanding up to leading order in powers of tjq, g 
and c, one obtains the same expressions as ((18(1 and 1(20(1 , 
except that the last term is missing. 

The assumption that quenched surface disorder is 
equivalent to annealed surface disorder has been used in 
replica mean field calculations^ and Monte Carlo simula- 
tions^, based on the argument that the chain "samples 
all significant surface patterns with a probability that 
asymptotically approaches that of an annealed system in 
the adiabatic limit" (cited from Ref. i). Here we recover 
this equivalence in the limit of infinite surface area S. 

We turn to the calculation of 

[(S(x(n + k)- x')V{z{n k))5{x{n) - x")F(z(n))]o. 
We assume the validity of the ground-state dominance 
approximation for the overall chain and for the side 
subchains. This gives 



[<5(x(n + k)- yL')V{z{n + k))S{x{ n) — x")V{z{n))]o 

J dn dY2 drg dn G(ri, ra; n)5(r2 - (^'))G(r2, rg; k)5{vs - (^'') )G(r3, r4; N - n ~ k) 



(32) 



1 



/ dvi dv2 dr^ dr^ G(ri, rs; ri)G(r2, rg; fc)G(r3, r4; N -n-k) 

bi 



2fc2/fco 



Gx(x',x";fc)G.(zo,^o;fc)e-l^«l\ ^ , , 



We consider first the homopolymer case. Inserting 
into lUHl), we obtain 



2fc2/fco 



N 1 + {2ki, - ko)zo 

J dxdx' .g(x-x') ^(x, zq, x', zq, |eo|)| ~ Pvq 



Here 



(33) 



g{Y,Y',p) = / dke-P^Giry-k) (34) 

~ / dke-P''G{r,r';k) 
Jo 

is the Laplace transform of the "complete" Green's func- 



at z = z' = zq (i. e., at the level where the poten- 
tial acts on the chain). The procedure of calculating 
Q{x, zq,x' , zq;\£o\) is briefly described in the Appendix 
IbI The result can be written in the following integral 
form 



^(x, zo,x',zo; |eo|) 



6 1 



dq e 



-iq(x— x') 



a2 (27r)2 
sinh (zo\/g^ + kfj 



(35) 



where q is the two-dimensional vector and q = |q|. In- 
serting ((35(1 into ((33() and changing the order of integra- 
tion over q and x yields the following expression for the 
free energy difference between the reference system and 
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N 



1 + {2kb - fco)zo 



dq 



6 sinh I zq \fq 



- /3uo 



(36) 



The second case of interest is the neutral heteropoly- 
mer on the neutral heterogeneous surface. The difference 
between the Eqs. (|18|l and H20(l is due to the monomer- 



monomer correlation function Q). Inserting expressions 
into (EOJ), we obtain 



/3AF ^ 2fcg/fco 
iV I + (2fc6 - /co)^o 



dxdx' .g(x - x') ^(x, zo,x', zojTp + |eo|) - /3wo 



This finally leads to 
/3AF 



2fc2/fco 



TV 



1 + (2fcf, - fco)zo 



dq 



6 sinh (zo 



(q^ + r2)3/2 (^(^) g2„a;(9) _ fc^ sinh(zo ^{q))) 



(3vo 



(37) 



(38) 



with iu{q) = ^q^ + eTp/a^ + fcg. Note that (jSHI) can 
formally be obtained from H38|l by setting Ap = 1 and 
Fp — > 0. This is, however, not a general rule. The cor- 
responding expression for a heteropolymer adsorbed on 
a homogeneous surface, Eq. (41) in Ref. l2^ is not fully 
contained in Eq. (|S5|> ^ 



E. Transition point 



is desorbed in both the original and the reference system. 
For the other, nontrivial, solution we again consider both 
cases separately. 

In the homopolymer case, the nontrivial solution for vq 
is given by setting the term in square brackets in Eq. (|36|l 
to zero. The transition point for the reference system cor- 
responds to the condition — zq^o = 1 and fc^ — *■ 0. Sub- 
stituting these conditions and introducing, for the sake 
of convenience, the dimensionless variables 



Having obtained an approximate expression for the dif- 
ference A_F between the free energies of the original and 
the reference systems, we adjust the auxiliary parameter 
vq such that the free energy of the reference system is the 
best approximation for that of the original one. Thus we 
require IAF{vl) = 0. For the Ai^ given by or P^ . 
this equation has one trivial solution fcb = or = 0, 
corresponding to the uninteresting case of a chain which 



f3vo 
a 



CTO 



a 



Ts = aTs, and zq = — (39) 

we obtain an equation for the transition point cto^o)*'^ in 
the random system; 



-l3vo = -6 A, ^2 / 



sinh (zoq) 



g2 



3/2 (g^og _ sinh(zog)/(2o(7)) 



(40) 



where (i5o^o)homo = l/(6-So) is the transition point for the 
homopolymer near an attractive surface. Eq. (|40|l de- 
termines the relation between a, Fg, and A^ that corre- 
spond to the adsorption-desorption transition in the orig- 
inal heterogeneous system. 



We turn to the case of a neutral heteropolymer on a 
neutral heterogeneous surface. Here, the nontrivial so- 
lution for vq, is obtained by setting the term in square 
brackets in Eq. (|38|1 to zero. We insert again the con- 
ditions for the transition point in the reference system. 
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kh = and fco zo = Ij ^md switch to the dimensionless dimensionless. Since (Tq = = 0, the equation for the 
variables according to Eq. (|39|l . The parameters Tp transition point then reads: 
and Ap associated with the RHP statistics are already 



dq 



6 sinh (zoi/g^ + Q^p) 



'A2 

p 



GFp e 



(41) 



IV. SUMMARY AND DISCUSSION 

We have used a reference system approach to calcu- 
late the adsorption/desorption transition of polymers at 
surfaces in systems with spatially correlated randomness. 
The resulting expressions become very elegant in the limit 
Zq ^ 1, i. e., if the range of the surface potential is short 
compared to the statistical segment length of the poly- 
mer. We shall now compile these results, contrast them 
with the corresponding result in our earlier paper and 
discuss them. 

We briefly recall the meaning of the different param- 
eters. ^0 denotes the mean polymer charge per segment 
and oq the mean surface charge per site; thus the mean 
affinity between polymer segments and surface sites is 
given by the parameter Q = ((Tq^o)- Randomness shifts 
the critical affinity parameter at the adsorption transi- 
tion Q'-*''-' to lower values, compared to the corresponding 
value (5ho2o of a homogeneous system. The shift depends 
on the variances Ap and of the polymer and surface 
disorder, and on the cluster parameters or inverse corre- 
lation lengths Tp and Tg- 

The case of random heteropolymers on homogeneous 
surfaces was already considered in our earlier paper—. 
The transition in the random system is shifted with re- 
spect to the homogeneous system by 



Q 



(tr) 




(42) 



The case of homopolymers on random surfaces is covered 
by Eq. In the hmit zq ^ 1, the integral in can 

be calculated explicitly and the solution takes the form 

6 



q£2 



^2 „ 



(43) 



Heterogeneous polymers on heterogeneous surfaces were 
discussed in the special case of zero mean affinity. The 
result, given by Eq. (|41|l . simplifies considerably in the 
limit zq 1. The heteropolymer is adsorbed, if 



6. 



A, A 



> 



const. 



(44) 



This defines a phase boundary in the parameter space of 
variances Ap and Ag and cluster parameters Tp and Tg. 



In order to complete this compilation of results, we 
also give the general expression for the shift of the ad- 
sorption transition in systems of heterogenous on het- 
erogeneous surfaces (not necessarily neutral). It con- 
tains the sum of the three contributions H42|) . (|43|l . 1)44(1 . 
but has one extra term from the expansion of = 

J^kLoilo^ dni g(ni,n2) c(n2 - n-i))'' in CHI)- The cal- 
culation is absolutely analogous to the one presented in 
Sec. (jnH) and yields 



-Alir^J^-A.^ll (45) 



^a/a2c + o(a/,a4) 



where the coefRcient C of the higher order term As A^ 
is given by 



C ^ 




with the Meijer G function G. The calculation of C is 
shortly sketched in the Appendix lO Note that Eq. (|45() is 
an implicit equation for Q^*^\ because the affinity Q = 
(o^Co) depends both on the mean charge per segment 
of the polymer and the mean charge oq per site on the 
surface. 

The validity of the approximation zq ^ 1 is 
demonstrated in Fig. ^ which shows some adsorption- 
desorption transition curves for homopolymers adsorbed 
on heterogeneous surfaces at different values of zqj as 
calculated from Eq. H4U|) . As a function of the correla- 
tion length l/Fs, the curves become straight lines with 
increasing 1/F with the same slope as the asymptotic 
curve, Eq. H43|l . The slope of the curves deviates from 
the asymptotic value only at small correlation lengths. 
The simplified expression H43|l approximates the full solu- 
tion l(in|l reasonably for zq < 1- Thus the approximation 
is justified if the range of the potential is shorter than 
the statistical segment length. At first sight, this con- 
dition seems somewhat arbitrary, given that the value of 
the statistical segment length a depends on the particular 
choice of the "segment" unit. One can always increase a 
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2 4 6 8 10 



1/r 



FIG. 1: The shift (Q"' - (5L''omo)/As of the adsorption- 
desorption transition point for a random surface with expo- 
nentially decaying sequence correlations compared to that of a 
homogeneous attracting surface, vs. correlation length l/Fs, 
according to Eq. 1401 1. The different curves show the results 
for different values of £o as indicated. The bold line corre- 
sponds to £o = 0, Eq. 1431 . 



by choosing a more coarse grained polymer description, 
i. e., combining several segments to one new effective seg- 
ment. The real length scale set by the non-adsorbed poly- 
mer is its radius of gyration. Applying the approximation 
zo 1 in the calculation of the adsorption transition is 
thus justified whenever the range of the potential is much 
shorter than the radius of gyration of the polymer. 

Our results can be summarized as follows: The pres- 
ence of correlations always enhances adsorption. Both in- 
creasing the correlation lengths, i. e., decreasing in Tp or 
Fg, and increasing the strength of the correlations (given 
by Ap, As) results in a shift of the adsorption transition 
towards lower affinities. 

The results obtained for homopolymers on random sur- 
faces ()43ll and those for heteropolymers on homogeneous 
surfaces (|42|) are very similar. The only difference are 
the exponents of F,, and Fp, which reflect the different 
"dimensions" of a planar surface and a linear polymer 
chain. 

The final result (|44|l for the adsorption of neutral het- 
eropolymers on neutral heterogeneous surfaces gives some 
particularly interesting information. First, it is worth 
noting that in Eq. H44() . both F^ and Fp enter in the de- 
nominator with their "characteristic" exponents (1 and 
1/2, respectively). The other notable feature is that As 
and Ap enter in (|44|l as the product of their squares 

— ^2 

As Ap - a similar dependence on the variances has been 
observed by Srebnik et al. — for the case of short-range 
correlated RHPs at a random surface. This means that 
if the surface (heteropolymer) contains attractive, re- 



pelling, and neutral sites and is kept neutral on average, 
it is preferable to have strongly repelling and strongly 
adsorbing sites rather than weakly repelling and weakly 
adsorbing - "almost neutral" - sites. 

We note another interesting property of the considered 
model. In Sec. lIIIDl we have already discussed the equiv- 
alence of quenched and annealed surface disorder. The 
argument was based on the fact that all terms related to 
different replicas in Eqs. 118|) and (|2l)|l contain an extra 
factor 1/ S, which vanishes in the limit of infinite surface 
area S. Technically, the factor \/S reflects the loss of 
translational freedom of a replica, if it is restricted to 
stay close to another replica by a coupling through the 
correlation function (7(x — x'). The argument is thus only 
valid for surface disorder. In general, quenched and an- 
nealed disorder on the polymer are not equivalent. 

At the transition point, however, terms related to dif- 
ferent replicas always vanish, because every indepen- 
dent repHca also contributes a factor 2fc^/fco/(l + (2fcb — 
ka)zo) ~ 2fcf,. In other words, our model calculation 
predicts that the adsorption threshold is the same for 
quenched and annealed surface disorder. We believe that 
this is a general feature of the adsorption transition in 
disordered systems. Joanny finds the same equivalence 
when studying the adsorption of polyampholytes on sur- 
faces'^ within the Hartree approximation. He claims that 
it has been conjectured for similar problems by Nieuwen- 
heuizen2&. Indeed, polymers and surfaces have only very 
few contacts at the adsorption transition. It is conceiv- 
able that the contacts can be arranged in an optimal way 
regardless of the particular type of disorder. 

If this is true, it has obvious implications for the phe- 
nomenon of molecular recognition: It means that studies 
of the adsorption transition in systems with quenched dis- 
order, averaged with general, unspecific site distribution 
functions such as Eqs. © or ((S)), can only provide very 
limited information about specific recognition processes. 
The adsorption threshold of an annealed polymer, which 
has had the opportunity to adapt its sequence to the sur- 
face, is absolutely equivalent to that of a completely un- 
adapted (quenched) polymer. 

We can also use our results to discuss the issue of the 
interplay between cluster sizes on the surface and the ho- 
mopolymer. Fig. |5] shows an example of an adsorption 
phase diagram, obtained from Eq. H44|l . for neutral het- 
eropolymers on neutral surfaces as a function of the clus- 
ter parameters Fp and Fs. The phase boundary does not 
reflect particular affinities between certain cluster param- 
eters Fp on the polymer and Fs on the surface. Regardless 
of the value of Fp, the adsorption is always highest at low 

Fs and vice versa. We do not observe any complementar- 
ity effects. 

To summarize, we have been able to derive simple and 
elegant equations for the adsorption threshold in systems 
of polymers on surfaces with correlated quenched and an- 
nealed disorder. These may be useful for future studies. 
We have discussed our results in the context of molecular 
recognition and found that the complex process of pat- 
tern recognition between molecules and surfaces close to 
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Adsorbed 



Desorbed 



r„ 



FIG. 2: Adsorption transition of neutral heteropolymers on 
neutral surfaces as a function of Fp and Fs. The parameters 



are As 



1 and gJi'ri 



1. 



the adsorption threshold cannot be understood in terms 
of a simple model with quenched or annealed disorder. 
We conclude that recognition either takes place at higher 
affinities, beyond the adsorption threshold, as suggested 
by Chakraborty et al. — , or has to be studied by more 
sophisticated models and methods. For example, the sit- 
uation might change if the probability distribution itself 
reflects the pattern on the surface. Such considerations 
shall be the subject of future study. 
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and and calculate 



/ = 



/I?{C(n)} exp(ea + \m eM^W - 0^'^^ - 0) 
exp 

/I?{C(n)}e-^«['="'-51« 



/PU(n)}e-^ 



(Al) 



The denominator corresponds to the normalization fac- 
tor in Eq. and the expression j 'D{^{n)} denotes 
functional integration in sequence space. To sim- 
plify the first exponential on the right hand side of 
Eq. (im . we define K = 1 - gc (Eq. (HII) and 
replace [c~^ ~ 9]~^ = cK^^ thus obtain the expres- 
sion exp + acK-^a/2 + ^K-'^a] . This cor- 
responds to the last three exponentials in Eq. H10|) . 

To calculate the integrals over sequences {$(«)} in 
Eq. ISU, 

we set n back to be a discrete variable, 
/ Vi^n)} = j;^ da • • • C d^N. This gives 



-SI? 



/P{eN}e^3?[^__ 
/2?{C(n)}e-^«c-^? 

1 



/ det(c-i) 
det(c-i -g) 
1 



(A2) 



^det[(c-i -g)c] ydctpq 
The final result, which enters Eq. H10|l . thus reads 



1 



exp 



VdetK 

In matrix notation, it is easy to sec that 



(A3) 



ln(dct(A')) = ln{det{T KT~^)) (A4) 
= ln([] A„) = ln(A„) = Tr(ln K) , 



APPENDIX A: PERFORMING THE SEQUENCE 
AVERAGE 



where T diagonalizes K, and A„ are the eigenvalues of K. 



We want to calculate surface disorder average in Eq 
lO. Let us define, for convenience. 



APPENDIX B: COMPLETE GREEN'S 
FUNCTION 



(j{n) ^ (3(Jo^V{za{n)) , ^{n) = 



and adopt the matrix notation 



N 

^G^ Ij^ dm dn2 e(ni) 



G{ni,n2) £,{n2) := / dn^{n)v{n). 

Jo 

The matrices c, c~^, and g are symmetric. In order to cal- 
culate the sequence average, we need to combine Eqs. dHl 



Our starting point is the equation H21() with the initial 
condition (|22|) and the boundary conditions f l23l) and 
1)24(1'): Laplace transforming with respect to n and taking 
into account the initial condition H22|) . we obtain 

-p g(x, z, x', z';p) + 5(x - x')(5(z - z') (Bl) 

= — g-V^CJ-/3wo S{z-zo)g{x,z,x.',z';p), 

where C/(x, 2, x', z';p) is the Laplace transform of 
G(x, z, x', z'; n) with respect to n with the boundary con- 
ditions 



g{x,z,x.',z';p)\ 
g(:ic,z,x.',z':p)\ 



X| .Z — 'OO 





0, 




(B2) 
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Firstly, we apply to this equation a Fourier transform 
with respect to x (x ^ q and Q ^ 



-p^(q,z,x',z';p) +e-i^(5(z- z') 
^2 ^ , .o.. OH 



—(3wq5(z — zo)^(q> z.x.', z';p) 

and secondly a Laplace transform 
respect to z (z — > s and - 



(B3) 



with 
.9)- 

I 



-Pg{ti,s,x.',z';p) + e '1'' 



6 



[(-|q|2 + ,s2)g(q,s,x',z';p) 



(B4) 
(B5) 



-.s ^(q, 0, x', z';p) - ^'(q, 0, x', z';p)] 
- Pwo e-^^'giq, zo,x',z';p) 

Here denotes d^/dz. This equation can be solved for 
.g(q,s,x',z';p), giving 



5(q, s,x',z';p) = 



2\ -z's -qx' 



^^(q,0,x',z';p) - (6/3wo/a^)^(q, zo,x^z» e"^"^ - {6/a^)e-' 

— |qp — dp/a^ 



(B6) 



Taking the inverse Laplace transform of ljB6|) with respect to s and the inverse Fourier transform with respect to q, 
we obtain for (x, zq , x' , zq ; p) : 



^(x, zo,x',zo;p) 



6 1 



dq e ^ ^ sinh {^Zq^ + 6^/0^^ 

(27^)^ V ^g2 + 6p/a2e"nV9'+6p/'^' - /co sinh (^zo^/^fT&p0) 



(B7) 



r 



APPENDIX C: CALCULATION OF THE 
HIGHER ORDER TERM IN EQ. HH 



this expansion corresponding to the order Af will have 
the form: 



Using the matrix notation of the Appendix ^ we ex- 
pand as 



(CI) 



fc=0 



exp 



alii / dnidn2dn3g{ni,n2)c{n3-n2)(3y^^ V{za{ns)) 
•^0 



(C2) 



This is equivalent to the following contribution to the 
and insert it into (|10|) . Then, the additional term from (3AF/N before taking the replica limit m 0: 



N 



f3^<^Uo dx'dx"g{x'-x") / dnidn2dn3.c{n^-n2) ^ y(z„(ni))(5(x„(7ii) - x')F(z^(n2))(5(x^(n2) - x")F(z,,(n3)). 

m 



r 



Here, " denotes the average in a reference system of m 
independent homogeneous replicated systems. First of 
all, we notice that only m terms corresponding to replicas 
with the same indices (i.e. a = 7 = 77) will contribute 
to A_F in the transition point. After taking the replica 
limit, we thus have 

(3^ alii j dx' dx" g(x' — x") J dni dn2 dn^^ 0(713 — n2)x 

[y(z(ni))<5(x(7ii) - x')F(z(n2))(S(x(n2) - x")F(z(n3))]o 

(C4) 

Note that the monomers ni and n2 are "coupled" in 



the expression (|C4|) via the site-site correlation function 
g, whereas the correlation between monomers 712 and 713 
is accounted for by the monomer-monomer correlation 
function c(?Z3 — 712). Therefore, one has to take into ac- 
count the order in the sequence of monomer numbers ni, 
7^2, and 713. There are 3!=6 possible permutations of rii, 
712, and 713 but in fact it is sufficient to consider the three 
possible arrangements schematically shown in Fig. 13 The 
other three are inverse of these ones. We calculate the 
average in Eq. (|C4p using the Green's function of the ref- 
erence system, like we did in Eq. H32|) , and then integrate 

it over J dx' dx" and J^dni dn2 dn^ . These calculations 
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n, < < n, < n, < n, < < 



(a) (b) (c) 

FIG. 3; Diagrams contributing to the expression IIC4t . See 
text for explanation. 



where G is the Meijer function. Taking into account the 
double counting (due to inverse sequences) and the pref- 
actor (To ^0 we arrive at the second term in the expression 
gSJ for the coefficient C in Eq. gsj . The first term in C 
comes, as it has been aheady noticed, from the expansion 



are straightforward and give the foUowing resuhs (in the 
dimensionless variables Ts and As): 



(a) 
(b) 

(c) 



A^A. 



-2 36 1 



2 36 1 



/6r„ 



2~2 36 22 

ApAs _2 ^22 y 



7rr„ 



6r, 



1 1 

1/2 3/2 



(C5) 
(C6) 

(C7) 
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